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We discuss the effects of surfaces on bulk plasma oscillations in a metallic film of 
thickness L. The leading corrections to the plasma frequencies due to the surfaces are 
proportional to L^ 1 . The frequencies of low-lying long wavelength modes are given 
by the dispersion relation of bulk modes, u = uj p [1 + a(q 2 /kp)], where the allowed 
g-values are = (nir/L), L = L — 2<4, and c4 is a complex length characteristic of 
the surface (reminiscent of but different from Feibelman's d±). An explicit expression 
for d;, is derived. Resonant excitation by an external field, E§e~ %ult is calculated. 
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The widespread interest in nanostructures (thin films, layer structures, quantum wires, 
quantum dots etc.) as well as in atomic and molecular clusters, in which at least one 
dimension is "mesoscopic" (10 A - 1000 A), has led us to re-examine the effects of finite size 
on charge collective modes (plasmons). In an infinite system, these modes have frequencies 
given by the bulk plasmon relation [|TJ 

u(q) = u p + a (J^J u p , (1) 

where u p is the classical plasma frequency, ui p = 4iie 2 n/m, with n the electron density, k F 
the Fermi wavevector, and a a (complex) number. The imaginary part of a describes the 
damping of bulk plasmons, which in an infinite system is due to excitations of multiple 
electron-hole pairs. In a finite system, we expect the following qualitative changes in the 
plasmon spectrum due to the presence of surfaces. First, only a discrete set of frequencies 
will be allowed. Second, there will be additional damping due to the decay of plasmons into 
(single) electron-hole pairs. In this Letter, we will show how both the real and imaginary 
parts of the allowed frequencies are completely determined by a single (complex) length. 
Specifically, we consider a thin slab with a positive jellium charge background, n + = n 
between z = and z = L and infinite in the x and y directions, and a neutralizing electron 
liquid described by the equilibrium density distribution n_(z) (See Fig. [I]). We shall discuss 
bulk-like collective modes of the form 

n(r,t) = n{z)e- tuj(n) \ (2) 

where the eigenfrequencies uj^ are near the classical bulk plasma frequency associated with 
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the density n, ui p = (47re 2 n/m) 

This problem has been previously addressed experimentally |2||| and theoretically [|J for 
modes of the form @ and of the more general form 

n(r,t) = n(z)e i( * x - ut \ (3) 

In the present Letter, limited to the case p = 0, we derive the following new result: the 
(complex) eigenfrequencies u^ n \ near tu p , are determined by the bulk plasmon dispersion 
Eq. (HD for small q, and by an effective, complex, slab thickness 

L — L — 2d b , (4) 

where dj,, precisely defined below, is an effective, complex, surface thickness. The thickness 
db is reminiscent of the parameter d± introduced by Lang and Kohn for uj = and 
generalized for arbitrary uj by Feibelman but is not the same quantity. The allowed 
(complex) interior wavenumbers, q^ n \ are given by 

q in) L = nn (n = l,2,...) (5) 

and the corresponding (complex) eigenfrequencies are given by Eq. (JJ): 

w (n) =oo(q in) ). (6) 



2 



We now derive the above-stated results. We take the equilibrium electron density, n-(z), 
to vanish for z < — c and z > L+c, and introduce an auxiliary thickness a, large compared 
to the surface thickness but much smaller than L. We then distinguish two surface regions, 
—c<z<a and L — a < z < L + c, and the bulk region a < z < L — a, in which n_(z) = n 
(see Fig. §. 

We consider a collective mode of frequency u, with self-consistent electric field, E(z), 
parallel to z. This field induces a current density, 

j(z) = dz'a{z,z')E{z'), (7) 

where a(z, z') = J dx' dy' o zz {x — x',y — y';z,z'); a zz is the zz-component of the non-local 
conductivity tensor. (For ease of notation, the time dependence, e~ lu)t , and the functional 
dependence of fields, currents etc. on uo will generally not be explicitly written out.) The 
quantity <j(z, z') is short-ranged in the difference variable (z — z'); also, when both z and z' 
are in the bulk region, 

a(z, z') = cr(z — z'), a < z, z < L — a (8) 

equal to the bulk conductivity corresponding to n and u>. Of course <r(z, z') vanishes when 
either z or z' are outside the slab (— c, L + c), as do E(z) (due to charge neutrality) and 

Denoting the induced charge density by n(z), the continuity equation and Gauss' law 

are 

iuen{z) + = (9) 



and 

dE(z) 



dz 



-Anen(z) (10) 



(where e has been taken as positive). Since at z = — c both j(z) and E(z) vanish, (|[) and 
(|10|) imply that, for all z, 

j(z) = 1 £e(z). (11) 

We can eliminate j(z) from (|7|) and (^1|) which results in the integral equation 

t £E(z) = J^ C dz'a(z,z')E(z'). (12) 

We are looking for solutions, which, in the bulk, behave as e ±iqz , with q -C kp. In this 
limit of long wavelengths in the bulk, we can, with Eq. (|5p in Eq. ( |T2"D obtain the frequencies 
of these solutions by expanding E(z') about z when z and z' are in the bulk. The result is 

^ = a (uj) - ^a 2 {u)q 2 , (13) 
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where 
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/1UJ 
dz'a(z-z';u) = ^, (14) 

a 2 (u) = J dz'{z' - z) 2 a(z- z';u). (15) 

Equation (|I3| ) is the bulk plasmon dispersion relation of an infinite system of conductivity 
a(z — z'). In the limit L 3> a, the frequencies of the lowest-lying eigenmodes of Eq. ([12]) are 
very near u p and the shape of E(z) in the surface region tends to the limiting shape 

E(z) -> E(z) = E(z; u p ) (16) 

shown schematically in Fig. EJ We arbitrarily take E(a) = 1 and write, near z = a 

£ w ^ (17) 

where ofo is the, as yet undetermined, point where E(z) extrapolates to 0. When z is in the 
bulk the solution which joins (0) near z = a is 

E(z) = , 1 ,. swq(z-d b ). (18) 

Thus, df, plays the role of a (complex) scattering length. 

The length db can be simply expressed in terms of the limiting surface density distribu- 
tion, n(z), associated with collective modes with q — > and uo — > u; p . From Eq. ([H]) and 
the form of E(z) in the surface region we see that n(z) has the general form shown in Fig. 
El. Thus, for z near a, we have 



E(z) 



dz' n(z) = / dz' 
Aire) J~c Jo 



d ,dn{z') 
-{zn{z)) - z 



dz' dz' 



- , , rfn(z') 

zri-r — I dz z — 

dz 1 



oo 



= zl dz'^l- Tdz'z'^1. (19) 
J-c dz' J-c dz' 

The upper limits oo are to be understood as ^> kp 1 but <C L. The linear form of E(z) near 
a, Eq. (|T9|) , extrapolates to zero when z equals 

4 = U /* ■ (20) 

/•oo » rfn(z) 

Thus, dft is the center of mass of the function dn(z)/dz. 

The eigenmodes are either even or odd about L/2 so that, by Eq. fll8|), the allowed 
values of g are 

g (n) L = mr, (n=l,2,...) (21) 
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where L is give by Eq. (|j). The corresponding frequencies are given by Eq. ([!]). Using the 
notation d b = d' b + id' b ' etc, we have 



k 



H a " +4 4)(i) 2 - « 23 » 



where L' = L — 2d' b . We observe that the real parts Eq. (|2"2"| ) of the frequencies correspond to 
bulk plasmon modes with the boundary conditions E(d' b ) = E(L — d' b ) = 0. The imaginary 
parts, (P^D, have contributions from the bulk (oc a") and from the surface (oc (d b /L )). 
We have solved the integral equation ([12]) in the discrete form 
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—E m = b) y a mn E n} (24) 
47r n=0 

where b = L/N is the interval length. We have chosen L = 35.4 A, N = 248, b = 0.14 A, 
and c = 0. For the sake of illustration, we have chosen the following values for a mn . In the 
"bulk" , 4 < m < 125, we have chosen 

■ ( \ iuJ p u p r • i 1 • 

Om, m = in(to) = - r - - , 2 ,o «2 - «*i , cr mm+2 = cr m m _ 2 = 

Attujo Snkpb 6 2 

2,U!p . 1 . , , 

and symmetrically for m > 125, with fcp = 1.75 x 10 8 cm -1 , hu p = 15 eV, and a = 
0.57 + i0.025 in Eq. ([!]) (appropriate [[[J for Al in the limit of q/kp <C 1). The functions 
niyj) and z/ are chosen so that the average and second moment of a mn are aoiuj) and a2{w p ). 
In the "surface", < m < 3, we have chosen 



O"0n 
&m,m 


= <7 n0 = 


+ e; 
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m,m+l 




m 

7' 


>V + e) 


m,m+2 




m 
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>V + e) 


m,m+3 




m 
16 


(?V + e 



cr 7l 

C"m,m+3 = 0"m+3,m = (^ + e) (26) 

and symmetrically near the other surfaces. The small real part e = u p a" / (6nk F b 3 ) was added 
to the conductivity in the surface regions as a simple model of damping in the surfaces due 
to decay into electron- hole pairs. Apart from this small real part, equation (|26"D is a simple 
interpolation between o" n = cr n o = and the bulk conductivity. The results for the lowest 
three modes are plotted in Figs. |] and f|, and confirm our conclusions, Eqs. (^[), @, ( |22| ) 
and (0). 

How do our results relate to charge density waves in classical Maxwell theory? "Classical" 
here means a local a(z, z') 
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Odass^, Z 1 ) = (T 5(z - z') = ~r^-5{z - z'), 

4:7100 



(27) 



so that (|T2l) becomes 



fra- 



cas) 



The boundary conditions |§ are, from charge neutrality, E(0) = E(L) = 0. Thus, classically 
the general solution is 00 = u p and E(z) = F(z), an arbitrary function of z satisfying the 
boundary conditions E(0) = E(L) = 0. In particular, the functions 



FJz) 



®n(QdLz) QciL = ^ n = 0, 1, . . . 



» 



77,7T 



(n) 



tOp, all n, 



(29) 
(30) 



are solutions. Thus, our non-local and translationally non-invariant a(z, z') leads to: (1) 
replacement of L by the complex L = L — 2<4; (2) a /mite dispersion of u^ 1 '; and, if cr(z, z') 
is known in detail in the surface region, the detailed field and density distributions in the 
surface region. 

We now turn to resonant response to a uniform external field, E e~ tuJt . Eqs. (|7|), ([|) 
and fllPP remain unchanged, but in Eq. (|TTD E(z) is replaced by E(z) — E and Eq. flTjp 
becomes 



7'/,; rL+c 

— [E(z)-E } = J dz'o-{z,z'-u)E{z'). 



(31) 



Near w = w p we can write this as 







+ 77 




4tt v ; y 



,r(z,z') 



47T 



-E'O) 



(32) 



where (7(2;, z') = cr(z, z'; cj p ), r(z, 2/) = [dcr(z, z'; u)/duj} w=uj and r] = uj — uj p . On the right- 
hand side, 77 can be neglected. 

Setting Eq = 0, for a moment, gives us the equation for the slightly complex normal 
modes, E n (z), and eigenvalues, rj n = 7]' n + ir]^. The kernels a(z, z')/i, t(z, z')/i are symmetric 
and (slightly) complex. The E n (z) satisfy the orthonormality relations 



/.L+e rL+c 
J dzdz'E m (z) 

We now substitute the Ansatz 



47T 



8(z - z') 



EJz') 



E(z) = J2A n E n (z) 



(33) 



(34) 



into (||), giving 



Y,A n (r ] -r ]n )E n (z) = ^E . 

47T 



(35) 
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Using the orthonormality relation Eq. ([33]) leads to 



A r . 



uj p E 



An(r] - r] n ) 



/L+c rL+c 
dz / dz'E n (z' 
-c J-c 



U{z - z>) - T -^l 
An i 



(36) 



Only even modes, n = 2,4... are excited. A well defined resonance occurs for n ~ rj' n if 
Vn ^ l^n — ^±2 1- I n this case as sm gl e term in ( |34"D is dominant, giving 



-En 



[(^ - V'n 



An 



(V ~ Vn) 2 + V'n 



(37) 



In order to resolve the discrete modes experimentally, the spacing between the modes must 
be larger than the width. From Eqs. ( p2|) and ( p3| ) we find that the spacing and width are 
given by 



(v'n ~ Vn-2) 



n 
Vn 



a 



n 



X 



(An - A) 




(38) 
(39) 



Thus, as n increases, the width gains on the spacing. For a film of width V = 30 A, 
d' b = d'l = 1 A, and a' = 0.57, a" = 0.025 (Al), the highest well defined resonance is n ~ 20. 

We plan to address the considerably more complex theory of modes with finite wavevector 
p (Eq. H) in the near future. 
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FIGURES 



FIG. 1. Positive background density, n + (z) (dashed line), and equilibrium electron density, 
n—(z) (dotted line). 

FIG. 2. The limiting (u) — > uj p ) form, E(z), and mode density n(z) in the surface regions. The 
mode density n(z) approaches the constant value no in the bulk. 

FIG. 3. (a) The real part of the electric field of the three lowest bulk-like eigenmodes near the 
surface. In the surface, the fields all approach the limiting form E{z). Insert The real parts of the 
three lowest bulk-like eigenmodes for < z < L/2. In the bulk region, the fields are sinusoidal. 
The vertical dashed line marks the position of a. 

FIG. 4. (a) Real and (b) imaginary parts of the plasmon eigenfrequencies. The straight lines 
are the dispersions Eqs. ( ^2| ) and ( |23| ) with V = 34.8 A and df," = 6 x 1CT 2 A obtained from Eq. 
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